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Graphical Solution for High- Angle Fire, 

By Prof. A. N". Whitehead, F.K.S., Imperial College of Science and 

Technology. 

(Eeceived November 20, 1917.) 

The point of the following methods is the reduction of the calculation of 
the path of a shell in high-angle fire to processes which can be either 
graphically constructed or calculated, once and for all, the results applying to 
all trajectories and shells of any ballistic constant. 

It is assumed that the law of resistance H(v) is largely empirical, and is not 
to be reduced to any neat mathematical formula. Furthermore the variation 
of the density of the atmosphere is taken account of, and its dependence on 
the height is assumed to be given by some empirical formula. 

Let W be the weight of the projectile, v its velocity, -v^ the angle between 
its direction of motion and the horizontal, viz., 0^\/r<90 o for the shell 
ascending, and — 9 0°^-v^<0 for the shell descending. The resistance of the 
air is (W/Cg) x K(V), C being the ballistic constant of the shell, corrected for 
the state of the atmosphere. Thus C will vary slowly as the height (y) 
changes. 




Eesolving tangentially and normally, the equations become 

Cg^-K^-C^sint 1 
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Put 



tanh z — sin i|r. 



(2) 



Then ty will be called the " angular elevation " at the time t, and z the 
" hyperbolic elevation." 

vol. xciv. — a. 2 B 
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Eliminate t and ^ between equations (1) and (2), obtaining 

^5f&^ - (Cg)- 1 x E (v) + tanh z. (3) 

ctz 

Let p be the density of the atmosphere at the position of the shell. Then 

we can put 

= B/p, (4) 

where B is a constant depending on the special projectile and on the 
standard projectile and standard state of the atmosphere. 
We assume that p is some function of y. Put 

P=f(y)- (5) 

Thus (3) can be written 

d log* v 



'r>e 



dz 
Also it is easy to prove that 



= B" 1 xf(y) x g- 1 E (v) + tanh z. (6) 



-^ = —a ^tanhz. (7) 

dz 



These equations determine v and y as functions of z, assuming given initial 
conditions. 

We proceed to explain a graphical method by which, given any simultaneous 
set (v, y, z) and also given the increment dz to z, the corresponding set 
(v + dv, y + dy, z + dz) can be determined. The point of the method is that it 
proceeds by the use of five curves, the same for all projectiles and all 
trajectories. These curves can, therefore, be reproduced once for all on 
suitably arranged sheets of the proper dimensions. The actual construction 
required to find a series of values of (v, y, z) proceeding by small given 
differences dz, is such as to make a very moderate demand on a skilled 
draughtsman. 

The five curves are 

7] = tanh z, 7] = g" 1 E (v), 

V = log e ^, v = ff" 1 ^ V =f(y)- 

When z (or yfr) is positive, the shell is ascending ; when z is negative, the shell 
is descending. It may, perhaps, be convenient for the draughtsman to treat 
the latter case separately. 

Then putting z = —z, equations (6) and (7) become 

rflo f«fo>. = tmhz f "B' 1 xf(y)xg" 1 B(v), (6') 

dz 



dy 



dz 



_ = —g 1 v 2 tanh z\ (7') 
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Thus the graphical modification for the descending shell, if treated separately, 
is of the slightest description, and need not be further alluded to. 

In Diagram I we use equation (6) and proceed from v, z, y and z-hdz to 
v + dv. With as centre the plane is divided into four quadrants, and the 
curves in each quadrant are referred to the bounding lines of that quadrant 
considered as positive axes for that quadrant. 



(J,uac6ia*%£ f/) 



Quackant (2) 



[f-irffri 




CLSU4 Of V 



(J^jzc&tant fa) 



Diagram 1. 

In Quadrant (1) the curve r) = tanh z is drawn, also ON" = z, Q.N = tanh z, 
01 = 1, ON' = z + dz (dz being negative). 

The point M on OY is determined from the second quadrant (as explained 
below), so that 

OM = B- 1 x/( 2 /)xf 1 EW. 

Then QL, in Quadrant (1), is parallel to NM. Thus 

OL = B~ 1 x/(y)x£*~ 1 K('y) + tan]i& 

In Quadrant (2) the curves 

V = 9" 1 R 0) and v = f(y) 

are so drawn that the horizontal axis in this quadrant acts both as axis of v 
and axis of y. Then, on the horizontal axis 

ON! = v, ON a = y, OP 2 ' = P 2 ¥ 2 =/(?/), OB = B. 

2 b 2 
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On the vertical axis, OMi = PilSTi, and P 2 'M is parallel to BMi. 

In Quadrant (3) the curve r\ = log e v is drawn, so that the horizontal axis 
is the axis of v, as in Quadrant (2) for one of the curves. Thus KNi = 
log e ^ = OK. 

In Quadrant (4), S is the point (z, tj), where tj = log e v. 

Then SS' is drawn parallel to LI, to meet the vertical through W, where 
OX 7 = z + dz, and dz is assumed to be negative, as will be the case. Thus S', 
of equation (6), corresponds to S, when z becomes z-\-dz. Then, working 
backwards, B/ is determined in Quadrant (3), and thence N"/, so that 

ON"/ = v + dv. 

Hence v + dv 9 corresponding to z -}- dz, is determined. 

Diagram II gives the graphical construction to determine y + dy in terms 
of z + dz, starting from the given simultaneous set (v, z, y). The quadrants 
are used separately as in Diagram I, but only two quadrants are now 
required, which (on the same plan as in Diagram I) are numbered Quadrant 
(1) and Quadrant (4) respectively. 




(J^uzdaa-ri^ // / 



Z-/ 



.nxtowk'Z 
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Chua£foa*z& f^-J 



"Diagram II. 



In Quadrant (1) the curves 

7] = tanh^ and tj = g~~ lr v 2 

are so drawn that the horizontal axis in this quadrant acts as the axis of z 
for the former curve, and as the axis of v for the latter curve. 
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Then on the horizontal axis 

ON x =:v, OPi 7 = P!Ni, 01 = 1, and ON = z, ON' = s + dfe, 

(dz negative in figure). 

On the vertical axis 

OQ' = QN = tanh'2, 

and L is determined by drawing P/L parallel to IQ'. 
On the horizontal axis produced backwards 

or == 1. 

In Quadrant (4), Y is the point (z, y\ and YY' is drawn parallel to PL, to 
meet the vertical through W in the point Y', which is (z + dz, y + dy). 

Thus v + dv and y-\-dy have now been determined in terms of v, z, y and 
dz. The procedure can be repeated, starting from (v + dv } z+dz, y-\-dy), and 
so on indefinitely. 

It is evident that it will not be necessary for a skilled draughtsman to 
draw many (or, indeed, any) of the straight lines which are here drawn for 
the purposes of explanation. 

It will be noted that the effect of height on air-resistance is taken into 
account in each step of the construction, and there is no smoothing out into 
zones. The effect of different states of the atmosphere on different days 
might be allowed for by choosing a sufficient number of such typical states 
and issuing different sets of skeleton sheets, each set corresponding to the 
form of the curve which belongs to that state of the atmosphere for which 
that set of sheets is to apply. 

Finally we note the form taken by the equation giving the dissipation of 
energy into the atmosphere, namely, 

d 1 

= (fy)-^K(tO J 

By this graphical method v and y have been tabulated in terms of z. The 
remaining elements of the path are given by 

1 f s 1 F 

t= I vdz, x = — iPsechzdz. (9) 

These quadratures can be effected either graphically as above, or by any 
other suitable method of approximation. Also y is given by equation (7) 
above. 

If a closer approximation is desired by an analytical method, we can 
proceed by using the graphical method as a rapid first approximation. 

It must be remembered that the whole point of any procedure, analytical 
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or graphical, is to reduce the solution as far as possible to processes which 
can be tabulated (or drawn) once and for all, suitable for all projectiles and 
all paths. 

Write u for the horizontal velocity, so that %i = v sec y{r. Kesolving 
horizontally and normally 

(I'll 

C — = — E(^ Sec i|r) COS "i|r, 
(X/v 

u ■—£- = — g cos"* y. 

(X/Xi 

Then [cf. equations (4) and (5) ] 

Bgu" 1 — =/(y)xE (u sec yfr) x sec ^. (10) 

For a certain range of relatively low values of v 

R (v) a v 2 , approximately ; 

for a certain range of high values of v 

R (v) oc ^ 3 , approximately ; 

for the intermediate range, no one value of n will give approximately 
R (v) oc v n . But, by dividing the whole range of important velocities for 
gun fire into a sufficient number of ranges, we can for each range find a 
value of n 9 say, n pq , which to sufficient accuracy allows us to write 

(Y p to Y q ), R (v) = fiv n p% 

These values of fi and n M can be found from the empirically constructed 
tables for R(v) by the method of least squares. Thus, if v = log e/ c6, and N"^ 
be the number of tabulations within the range Y p to Y q . 

N pq V -f 92^ S l0g e V = 2 l0g e R (tf), 
V 2 log e fl + ft^ £ {l0g e v} 2 = 2 {l0g e V . log e B(fl)}. 

The determination of n M for each range can be effected once and for all. 
With this value for n pq> 

R (u sec yfr) — R (%) {sec f%} 5 (V p to V ? ). (11) 

Then equation (10) becomes 

B ^ ST^ = 2 /(2/)(seof )'^ +1 ^f. (12) 

The integration on the left-hand side of (12) can be extracted from tables for 

C u du 
the function g -^y-y which is already tabulated in the Gunnery text-books. 
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Thus the only operation special to the path is the calculation of the right- 
hand side of (12). 

Now y is determined as a function of \{r by the preceding graphical method. 
Accordingly within each range ( "V^ to Y q ) or portion of such range, 

|/(2/)(seci|r)^ +1 ^ 

can be calculated by any convenient method of quadrature. Thus u (and v) 
are determined as functions of yjr, and all the elements of the path (viz., 
(%, y and t) are determined in terms of yfr. 



The Lommel- Weber fl Function and its Application to the 
Problem of Electric Waves on a Thin Anchor Ring. 

By J. E. Airey, M.A., D.Sc. 
(Communicated by Prof. J. W. Nicholson, F.E.S. Eeceived December 7, 1917.) 

1 f 1 " 
The function Q n (x) defined by the integral - sin (xsiiKfa—nfydcj), or, 

2 C w l 2 
when n is an even integer, by - sin (x sin <j>) cos n6d6 t is closely related to 

IT Jo 

1 [* 

J» (#), which was first given by Bessel under the form - cos (x sin $—n$) dfa 

7T Jo 

The X2 function, although of less importance in its application than the 

J function, occurs in a number of problems in mathematical physics, more 

especially those relating to the interference and diffraction of light. 

2 C n l 2 
The function O (x) = - sin (x sin 6) dd> was employed by Lord 

77 " Jo 

Eayleigh* in the problem of the reaction of the air on a vibrating circular 
plate, and by H. Struvef in the theory of diffraction in telescopes. H. F. 
WeberJ applied the function X2 f (a?) and, in addition, found a number of 
interesting properties of X2 v (x), including the recurrence formula and the 
development in series of ascending and descending powers of the argument. 

* Lord Eayleigh, ' Theory of Sound,' vol. 2, p. 164, equation (5). 

t H. Struve, "Beitrag zur Theorie der Diffraction an Fernrohren," c Annalen der 
Physik und Chemie,' vol. 17, p. 1013. 

J H. F. Weber, "Die wahre Theorie der Fresnel'schen Interferenz-Erscheinungeu," 
' Vierteljahrsschrift der Naturforschenden Gesellschaft in Zurich,' vol. 24, p. 48. 



